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The first high-statistics cross sections for the reactions -fp — > K* + A and jp — > A"* + £° were 
measured using the CLAS detector at photon energies between threshold and 3.9 GeV at the Thomas 
Jefferson National Accelerator Facility. Differential cross sections are presented over the full range 
of the center-of-mass angles, #2^+ , and then fitted to Legendre polynomials to extract the total 
cross section. Results for the K* + A final state are compared with two different calculations in an 
isobar and a Regge model, respectively. Theoretical calculations significantly underestimate the 
K* + A total cross sections between 2.1 and 2.6 GeV, but are in better agreement with present data 
at higher photon energies. 

PACS numbers: 



I. INTRODUCTION 

One motivation for the study of K* photopro- 
duction is to investigate the role of the if q (800) meson 
(also called the k) through f-channel exchange. The k 
is expected to be in the same scalar meson nonet as 
the /o(500) meson (also called the a). Neither of these 
mesons have been directly observed because of their large 
widths, which are nearly as big as their respective masses. 
Such a large width is expected for scalar mesons, which 
have quantum numbers J = ++ . In many quark 
models, there is virtually no angular momentum bar- 
rier to prevent these mesons from falling apart into two 
mesons, such as a — > tttt or k — > Kit. Because the a and k 
mesons cannot be observed directly, indirect production 
mechanisms provide better evidence of their existence. 

The a meson is rather well established [l[ as a 7T7t res- 
onance, which is an important component of models of 
the nucleon-nucleon (N-N) interaction such as the Bonn 
potential The K meson, however, is less easily estab- 
lished due to its strange quark content. Data for hyperon- 
nucleon (Y-N) interactions are sparse and hence models 
have a range of parameter space that may or may not 
include k exchange. Perhaps the best current evidence 
for the k is from the decay angular distributions of the 
Z?-meson into Kit-it final states H. 
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Here, in photoproduction of the K* + , the k + enters 
into the ^-channel exchange diagrams [|J. k cannot con- 
tribute to kaon photoproduction because the photon can- 
not couple to the K-k vertex due to G-parity conserva- 
tion. Theoretical calculations have been done [|| showing 
the effect of the k on photoproduction of K* + and K*° 
final states. Several years ago, two reports of K*° pho- 
toproduction were published BM 

but only preliminary 

results on K* + were available [7]. 

We present the first results of K* + A and K* + YP pho- 
toproduction with high statistics. Together, the K* + and 
A"* photoproduction results could put significant con- 
straints on the role of the k meson in i-channel exchange. 
Here, for the first time, we make the ratio of total cross 
sections for the reactions 7p —5- K* + A and jp —> K*°"E + 
and compare with the same ratio calculated from a the- 
oretical model for large and small contributions from k 
exchange. Other evidence for the k comes from recently 
published data on the linear beam asymmetry in photo- 
production of the 7p — > K*°T, + reaction 0] which shows 
a significant positive value at forward K* angles that is 
the signature of k exchange 

A secondary motivation for this study is to understand 
if theoretical models using Regge trajectories plus known 
baryon resonances can explain the K* + photoproduction 
data. If not, then there may be higher-mass baryon res- 
onances that could couple strongly to K*Y decay. In a 
classic paper on the quark model, Capstick and Roberts 
calculated many nucleon resonances that were pre- 
dicted, but not observed in existing partial wave analy- 
ses of pion- nucleon scattering. They also observed that 
some of the higher-mass resonances may couple weakly 
to pion decay channels and more strongly to KY and 
K*Y decays. Indeed, studies of KY photoproduction 
[Tol ] have shown that hadronic model calculations cannot 
explain the data without the addition of a new nucleon 
resonance near 1.9 GeV. We can look for other "missing 
resonance states at higher mass, such as those identified 
in the Bonn-Gatchina analysis [11| largely through pre- 
cise hyperon photoproduction data from CLAS, by com- 
paring K* photoproduction data to model calculations. 
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This paper is organized into the following sections. 
First, the experiment is described. Next, the data anal- 
ysis is presented in some detail. Then we compare the 
results with theoretical calculations. Finally, we discuss 
the significance of the comparison and provide some con- 
clusions. 



TABLE I: Some physical properties of the K* + , A and E° 
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A 



Mass(GeV) 0.89166 1.11568 1.19264 

Decay products K°ir + , K + tt° pir~ , mv A-y 
Branching fraction 66.7%, 33.3% 63.9%, 35.8% 100% 



II. EXPERIMENTAL SETUP 



The data used in this analysis are from part of the glla 
experiment, which was taken from May 17 to July 29, 
2004, using the CEBAF Large Acceptance Spectrometer 
(CLAS) located in Hall-B at the Thomas Jefferson Na- 
tional Accelerator Facility (TJNAF) in Newport News, 
Virginia. Real photons were produced by bremsstrahlung 
from a 4.0186 GeV electron beam incident on a 1 x 10~ 4 
radiation length gold foil. The electron beam was deliv- 
ered by the Continuous Electron Beam Accelerator Facil- 
ity (CEBAF). The Hall-B Tagging System [H was used 
to determine the photon energies by measuring the ener- 
gies of the recoil electrons using a dipole magnetic field 
and a scintillator hodoscope. The associated photon en- 
ergies were then calculated by the difference between the 
incident electron energies and the recoil electron energies 
with an energy resolution of about 2-3 MeV. The Hall-B 
Tagging System tags photons in the range from 20% to 
95% of the incident electron energy. 

A liquid hydrogen target was used in the glla experi- 
ment. The target was contained in a cylindrical Kapton 
chamber of 2 cm radius and 40 cm length. The target 
density was determined by the temperature and pressure, 
which were monitored once per hour during the glla ex- 
periment running. 

The CLAS apparatus was used to detect particles gen- 
erated from the interaction of the incident photons with 
the target. The CLAS detector was able to track charged 
particles that have momenta larger than ^200 MeV, and 
the detection area covered polar angles from 8° to 142° 
and 80% of the azimuthal region. It was composed of 
several sub-systems, arranged with a six-fold azimuthal 
symmetry. A plastic scintillator Start Counter, placed 
just outside of the target, was used to measure the ver- 
tex time of particles in coincidence with the incoming 
photon. The Start Counter was made of 24 scintillator 
strips with a time resolution of ~ 350 ps [ljj]. The su- 
perconducting coils of the CLAS detector generated a 
toroidal magnetic field that bent the path of outgoing 
charged particles. Those particles traveled through 3 re- 
gions of drift chambers [l4[ that measured the curved 
paths to give the particle momenta with a typical resolu- 
tion of ~1.0%. For the glla experiment, the current in 
the superconducting coils was set at 1920 A, which gave 
a maximum magnetic field of ~ 1.8 T. The time of flight 
(TOF) system was located beyond the outermost drift 
chambers at a radius of ^4 m from the target and was 
used to measure the time and position of each charged 



particle that hit the TOF scintillators. The TOF infor- 
mation, along with the particle momentum, was used for 
the particle identification in the analysis. The time res- 
olution of the TOF system was about 80 ps to 160 ps, 
depending on the length of the scintillators [HI . A more 
detailed description of the CLAS detector is given in Ref. 
10]. 

The event trigger for the glla experiment required that 
at least two tracks were detected in different sectors of 
CLAS. Once the event satisfied this condition, it was 
written to tape for future analysis. The data acquisition 
system for the glla experiment was able to run at ~ 5 
kHz with a typical livetime of 90%. 



III. DATA ANALYSIS 

As one of the largest photoproduction datascts at 
CLAS, the glla experiment has ~20 billion triggers. The 
calibration of each CLAS sub-system followed the same 
procedures as described in Ref. [17[. Additional details 
can be found in Ref. fl8| . 



A. Channels of Interest 

Because the K* + is an unstable particle, it will quickly 
decay to Ktt (see Table [TJ) by the strong interaction. 
By applying energy and momentum conservation, the 
K* + momentum is reconstructed from its decay parti- 
cles, K°7r+. The K° is a mixture of 50% K s and 50% 
Kl, but only the Ks decay is detected by the CLAS de- 
tector. The K° is reconstructed from the Ks decay to 
7r + 7r~ with a decay branching fraction of 69.2%. The 
same branching fractions are reproduced by the Monte 
Carlo detector simulations (see section 3.6), and hence 
are implicit in the detector acceptance values. 

To summarize, we report on the differential and total 
cross sections of the photoproduction channel(s): 

7 p -> A-*+A(£°) (1) 

followed by 

K*+ -»• K°tt+ (2) 
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and 



Ks — > 7T + 7T 



(3) 



The K* + and Ks are reconstructed directly from their 
decay products, while the A and E° are reconstructed 
using the missing mass technique. 



B. Particle Identification 

The Time of Flight (TOF) difference method was used 
to identify events with three pions (two positive and one 
negative charge) in the final state. Explicitly, 



Atof = tof n 



tofc. 



(4) 



where tof mea is the measured TOF of the particle and 
tofcai is the calculated TOF with the measured momen- 
tum p and the mass of a pion. In more detail, 



tof„ 



ttoi 



t. 



(5) 



where t to f is the time when the particle hits the TOF 
scintillators and t s t is the time when the photon hits 
the target. This information is determined by the CLAS 
Start Counter. In comparison, tof ca i is given by: 



where 



Thus 



f / L 1 

tO] ca l = — • — 
C p 



tofc 




(6) 



(7) 



(8) 



where L is the path length from the target to the TOF 
scintillators, c is the speed of light, p is the particle's mo- 
mentum, and m is the mass of a pion. The pion candi- 
dates are required to have \Atof\ < 1.0 ns. Fig. [T] shows 
the TOF difference spectrum. The solid lines define the 
region of the cut, the small peaks on the both side of the 
cuts are due to photons coming from other beam bunches, 
showing evidence of the ~2 ns beam bunch structure of 
CEBAF. 



C. Photon Selection 

After applying the \Atof\ < 1.0 ns cut, particles that 
came from different RF beam buckets were removed nat- 
urally. Of the photons measured by the photon tagger, 
we want those that come within 1.0 ns of the particle 
vertex time, which are called "good" photons. However, 
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FIG. 1: The TOF difference spectrum for pions. The two 
straight lines show the cut limits for selecting pions in a time 
window of ±1.0 ns. 



there might still be more than one "good" photon in each 
event. To select the correct photon, all "good" photons 
were scanned to find the one that gave the three-pion 
missing mass closest to the known mass of the A(S ), 
where 



MM (tt + it-7t + ) = J (&, + m p - - (. 



(9) 

is the missing mass summed over all three pions in the 
event, while E 7 and p 7 are the energy and momentum 
vector of the photon. The two-pion missing mass is sim- 
ilarly defined. 



D. Cuts applied 

Several cuts were applied to the data to reduce the 
background and to remove events below threshold for 
the reaction of interest. In general, the strategy is to use 
geometric and kinematic constraints to eliminate back- 
grounds while ensuring that the signal remains robust. 
The efficiency of various cuts was tested with Monte 
Carlo simulations (see section UlI Gp . 

The geometric and kinematic constraints used here are 
listed below: 

• Fiducial cuts were applied to remove events that 
were detected in regions of the CLAS detector 
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where the calibration of the detector is not well 
understood. 

• A cut on the vertex position along the beam axis 
(the z-axis) to be within the target position was ap- 
plied. All pions were required to be generated from 
the same vertex position within the experimental 
position uncertainty. 

• The missing mass from the K° was required to sat- 
isfy the relation MM(tt + tt^) > 1.0 GeV to include 
all hyperon mass peaks, for pion pairs with an in- 
variant mass inside the K° mass window (see next 
section). Similarly, the missing mass from the K* + 
was required to be greater than the nucleon mass, 
MM(7r+7r-7r+) > 1.0 GeV. 

After this step, the K* + A and K* + YP reaction 
channels were treated differently, since different 
backgrounds are present for each final state. For 
instance, the large background from 

A" £*+(1385) (10) 

present for the K* + A reaction channel makes the 
the extraction of K* + A yields by simply fitting the 
A peak impossible. On the other hand, there are 
only very small portions of the £* + (1385) that con- 
tribute to the K* + YP background, which can be 
easily removed based on Monte Carlo studies (see 
following section). Thus we could fit directly the 
S° peak in the three-pion missing mass for K* + TP 
channel, whereas a different approach (given below) 
is necessary to extract the K* + A yield separately 
from background due to JsT E* + (1385) production. 
The following lists the extra cuts applied for each 
reaction channel. 

• For the K* + A analysis, a cut was placed on the 
A peak in the three-pion missing mass: 1.08 GeV 
< MM(tt + tt^tt+) < 1.15 GeV. This ensures that a 
A was present in the final state. 

• For the K* + T,° analysis, a cut was placed on the 
K* + peak of the three-pion invariant mass: 0.812 
GeV < A/(7r+7r-7r+) < 0.972 GeV. This ensures 
that a K* + was produced. 



E. Sideband Subtraction 

Because reactions other than K* photoproduction are 
present, background is still mixed in with the channels of 
interest. Fig. [5] shows the two-pion invariant mass plot 
after the first three cuts in the previous section, inte- 
grated over all photon energies. A clear peak centered 
near 0.497 GeV sits on top of a smooth background. 
The invariant mass is calculated using the momentum 
vector of one tt + in the event, along with the 7r~ mo- 
mentum. Since there are two 7r + s, both 7r + 7r~ pairs are 
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FIG. 2: The reconstructed two-pion invariant mass showing 
the Ks distribution. The vertical lines define the bands for 
the SSM, as explained in the text. 

tested, but typically only one combination will satisfy all 
kinematic constraints. To avoid double-counting, in rare 
cases where both ir + satisfy all constraints, this combina- 
toric background is removed, for both data analysis and 
Monte Carlo acceptances. 

To reduce the background, a Sideband Subtraction 
Method (SSM) was applied. The concept of the SSM is 
to assume that the background in the signal region can 
be approximated by a combination of the left and the 
right regions, which are adjacent to the signal region. In 
our analysis, the two-pion mass of the Ks is used as the 
criteria to select the signal and sideband regions. Fig. [2] 
shows the regions used in our analysis. The middle band 
is the signal region, centered at the mass of K° with a 
width of 0.03 GeV. The other two bands, with the same 
band sizes, are the combinatorial background. 

Fig. El shows the sideband subtraction applied to the 
reconstructed three-pion invariant mass and to the three- 
pion missing mass. The SSM reduces the background, 
giving cleaner signal peaks. 



F. Peak Fitting 

After applying the SSM to each M(n + ir + ir~) invari- 
ant mass plot, corresponding to different incident pho- 
ton energy and different K* + production angle ranges, 
the K* + peak becomes clearer, but it is still not free of 
background. The main contribution to the background 
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FIG. 3: Three-pion invariant mass (left) and the three-pion missing mass (right). The four plots in each group correspond to: 
a) the left band, b) the middle band before the SSM, c) the right band and d) the middle band after the SSM. The peak in 
the three-pion mass is the K* + and the peaks in the three-pion missing mass are the A and E°. All plots are integrated over 
all incident photon energies. 



comes from the reaction channel jp — > K°I]* + (1385), 
which passed through all the cuts. In addition, the 3- 
body phase space reaction ^p — > K°ir + A is also present, 
and will contribute to the background as well. 

In order to extract the correct K* + peak yield, instead 
of fitting the K* + peaks directly with a Breit-Wigner plus 
background functions, we applied a template fit. The 
precondition for this template fitting is that we assume 
there is negligible interference between the K* + A and 
i4f°£* + (1385) channels; in other words, we assume that 
the K*+A and if°£*+(1385) add incoherently. If we re- 
move all other sources of background, then the K* + mass 
plot should have background only from the £* + (1385) 
peak. Similarly, the background in the E* + (1385) plot 
comes only from events in the K* + peak. Because the 
three-body K°tt + A channel is also a possible background, 
we assume it will add incoherently as well in both mass 
projections. 

To justify these assumptions, we explored the effect 
of various levels of interference between these two final 
states in the simulations. The result is that the template 
fits correctly reproduced the generated events to within a 
5% uncertainty for assumptions of maximal constructive 
or destructive interference. 

Fig. [4] shows an example of the template fitting, where 
the solid dots with error bars are from the data, while 
the curve is from the fit, which contains contributions 
from both the K*+A, A'°E*+ and K°tt+A channels. The 



K* + peak is seen in the left plots and the E* + peak is 
seen in the right plots. The template shape for each 
contribution comes from the simulation for that channel, 
and the magnitude of each channel is a free parameter 
to optimize the fit, with the result for each component of 
the fit shown in the bottom plots of Fig. |4j Both mass 
projections of Fig. U arc fit simultaneously to minimize 
the overall \ 2 - 

For the K* + T,° reaction, the counts from the E° were 
extracted by using a Gaussian fit, then the yields were 
corrected bin by bin based on a Monte Carlo study of how 
much if I]*+(1385) leakage there is to K*+Y,° reaction 
channel. The correction was studied and found to be 
less than 0.1%, which was included in our cross section 
calculation. Fig. [5]shows an example of the fitting. There 
are two peaks in the three-pion missing mass plot, one 
corresponding to the A and the other to the E°. The 
fitting function used two Gaussians plus a second order 
polynomial, for the A peak, E° peak and background, 
respectively. 



G. Detector Acceptance 

A computational simulation package, the CLAS 
GEANT Simulation (GSIM), was used for the Monte 
Carlo modeling of the detector acceptance. GSIM is 
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FIG. 4: Example of the template fitting. Top: the solid dots 
are from the data, while the curve is from the fitting, which 
contains contributions from the K* + A, K°Ti* + and K°tt + A 
channels, shown individually by the two plots at the bottom 
in large red diagonal cross, forward green diagonal and small 
blue diagonal cross histograms, respectively. 
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FIG. 5: Example of two Gaussians plus a second order poly- 
nomial fit to the reconstructed A and E° missing mass peaks. 



based on the CERN GEANT simulation code with 
the CLAS detector geometry. Thirty million — > 
K* + A (S°) events were randomly generated, with all 
possible decay channels of the final state particles (K* + , 
A, E°. . . ). The Monte Carlo files were generated with 
a Bremsstrahlung photon energy distribution and a tun- 
able angular distribution that best fit the K* data. The 
energy bin size was 0.1 GeV and the total cross section 
was assumed constant across the bin. This assumption 
is reasonable based on the slowly- varying total cross sec- 
tions shown below. Because the simulations have a bet- 
ter resolution than the real CLAS data, the output from 
GSIM are put through a software program to smear the 
particle momentum, timing, etc. to better match the real 
data. 

An extensive study of the glla trigger [13] showed a 
small inefficiency for the experimental trigger. To ac- 
count for the trigger inefficiency, an empirical correction 
was mapped into the Monte Carlo. The trigger correc- 
tions applied here is the same as used for other CLAS 
analyses of this same dataset [TtJ • 

The detector acceptance is calculated by: 



Dmc 
Gmc 



(11) 



where e represents the detector acceptance, Dmc is the 
number of simulated events after processing and Gmc is 
the number of generated events. 

The same software used for the experimental data was 
applied directly to the Monte Carlo data. Simulated 
events are extracted by fitting each reconstructed K* + 
peak for a given photon energy and K* + production an- 
gle. In our analysis, a non-relativistic Breit-Wigner func- 
tion 



1-4 



non — rBW 



* = a.L 



1 



2tt (E — Eji) 2 + r 2 /4' 



(12) 



was used to extract the counts of the K* + peaks for 
K* + A channel. Here, T is the full width at half max- 
imum of the resonance peak, E is the scattering energy 
and Er is the center of the resonance. 



As described in section MI D| different methods were 
used for the K* + A and K *+£° channels due to the pres- 
ence of A* resonance contributions in the former. For the 
K* + YP channel, where there is no kinematic overlap from 
hyperon states, the counts under the E° peak were fitted 
directly using a Gaussian function. Fitting the three-pion 
missing mass of the E° has less uncertainty than fitting 
the K * + peak, since the E° peak is relatively narrow on 
top of a nearly flat background. This method was used 
for both simulated and experimental data. 
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IV. NORMALIZATION AND CROSS SECTION 
RESULTS 

The differential cross sections are calculated by the for- 
mula: 

d ° Z fl3 l 

dcas&™ + N targe fN gflux -£-Acos6™.f u > K > 



where 



da 



is the differential cross section in the 



K* + angle center-of-mass (CM) frame, Y is the exper- 
imental yield, N target is the area density of protons in 
the target, N g fi ux is the incident photon beam flux, e is 
the detector acceptance, A cos 6^% is the bin size in the 
K* + angle in the CM frame and fu is the DAQ live time 
for the experiment. 

The detector acceptance e and experimental yields Y 
for the K* + A and K* + YP reactions are described in the 
previous sections. 

For each incident photon beam energy range (AE = 
0.1 GeV), nine angular regions were measured, uniformly 
distributed between —1.0 < cosO^li < 1.0. Hence, 
A cos «c« is §. 

The livetime fn for the glla experiment was carefully 
studied as a function of beam intensity, and found to be 
0.82 ± 0.01 for this measurement [l9[. 

In our analysis, photon flux was extracted in photon 
energy steps of 0.05 GeV. In the final analysis, we used 
photon energy bins of 0.1 GeV, and the fluxes added 
appropriately. 

The proton density N target is calculated using the for- 
mula: 



N, 



target 



p-L-N A 
A ' 



(14) 



where p, L and A are the target density, target length 
and the atomic weight of hydrogen, respectively. Na is 
Avogadro's number. For the glla experiment, an un- 
polarized liquid hydrogen target was used. The target 
density p was measured using: 



p = a\T +a 2 P + a 3 , 



(15) 



where T, P are the target temperature and pressure 
(measured at the beginning of each CLAS run), while 
o-i, 0-2, &3 are the fitting parameters. The mean value of 
the target density p for the glla data was obtained by 
taking the average (l7| : 



1 



P = 



N r 



p r = 0.07177 g/c 



(16) 



where N run is the number of runs. Using the target 
length of 40 cm, this gives N target . 



V. RESULTS 

Fig. [5] shows the differential cross sections for the pho- 
toproduction reaction jp — > K* + A, where there are 22 



plots, for Ej bins ranging from 1.70 to 3.90 GeV. There 
are nine angular measurements in each plot, uniformly 
distributed in cos 6*^+ between -1.0 and 1.0. In general, 
the K* + A differential cross sections shows dominantly a 
i-channel behavior, with an increase at forward-angles. 
Similarly, Fig. [7] shows the differential cross sections for 
7 p — > K* + S° photoproduction over the same photon 
energy range. Comparison with theoretical calculations 
are given below in section IV Bl 

The differential cross sections can be decomposed into 
Legendre polynomials as [Icj : 



do 
dcosd 



O total 



{1 



N 

£< 

(=i 



(17) 



where o'total is the total cross section. By fitting the 
differential cross sections up to 4 th order Legendre poly- 
nomials 



(18) 



i=0 



the total cross section was extracted by integrating f(x) 
over cos from -1 to 1. Using the properties of the Legen- 
dre polynomials, after the integration, only the do term is 
left. Hence the total cross section is given by a total = 2-ao 
Fig. [5] shows the fitting for the ^p — > K* + A channel, 
and Fig. [7] shows the fits for the K* + Y,° final state. The 
fitting parameters ao through 0,4 for each channel are 
plotted versus the incident photon energy E^ in Fig. HI 
The extracted total cross sections are shown in Fig. [S] 
for the K* + A and K* + TP final states along with some 
theoretical curves explained below. The error bars show 
only the statistical uncertainty. 



A. Systematic Uncertainties 

Systematic uncertainties come from several sources: 
the applied cut parameters, the choice of fitting func- 
tions, the Monte Carlo used for the detector acceptance 
and so on. 

Systematic uncertainties were estimated for each cut 
by varying the cut intervals and then recalculating the 
differential cross sections. The changes to cut parame- 
ters were applied to both the experimental data and the 
simulated output. The relative difference between the 
new cross sections and the original cross sections was 
calculated bin by bin using: 



So = 



&old 



Void 



(19) 



and then the resulting So values were histogrammed. 
This histogram was fitted with a Gaussian function, and 
the width from the Gaussian fit was taken as the sys- 
tematic uncertainty for each variation. The cut intervals 
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FIG. 6: Fitting the differential cross sections for jp — > K* + A with 4th order Legendre polynomials. Incoming photon energies 
range from 1.7 to 3.9 GeV. 



were varied to both larger and smaller values, and we 
chose the larger of the systematic uncertainties calculated 
from each variation. 

Similar estimation were done for the detector accep- 
tance, by varying the inputs to the Monte Carlo. Also, 
different fitting functions and background shapes were 
used to determine the systematic uncertainties associated 
with the peak yields. The total systematical uncertainty 



is then given by 

Stotai = ^51 + 51 + 51 + .... (20) 

which assumes no correlated uncertainties. 

The total systematic uncertainty from all sources, 
added in quadrature, is shown in Table [Til where the 
other sources include the target length, density and so 
on. For the K * + A final state the overall systematic un- 
certainty is 14% and for K* + T,° the systematic uncer- 
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FIG. 7: Fitting the differential cross section for 7p — > K* + T,° with 4th order Legendre polynomials. Incoming photon energies 
range from 1.8 to 3.8 GeV. 



tainty is 12%. 



B. Theoretical Calculations 

The models that are currently available for K* pho- 
toproduction are based on effective Lagrangians, which 
fall into two groups: isobar models and Reggeized meson 
exchange models. Isobar models evaluate tree-level Feyn- 



man diagrams, which include resonant and nonresonant 
exchanges of baryons and mesons. The reggeized mod- 
els, on the other hand, emphasize the i-channel meson 
exchange, which is expected to dominate the reaction at 
energies above the resonance region. The standard prop- 
agators in the Lagrangian are replaced by Regge prop- 
agators, which take into account an entire family of ex- 
changed particles with the same quantum numbers in- 
stead of just one meson exchange. In this section, the 
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FIG. 8: Legendre polynomial fitting parameters up to 4 th order plotted versus incident photon energy E-y for 7 p —¥ K* + K 
(left) and 7 p -* A"*+E° (right). 
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FIG. 9: Total cross sections of the reaction 7p — > K* + A (left) and 7p — ^ _K"* + E° (right). 



K* + A cross section results will be compared with calcu- 
lations from these two theoretical models. 

One model we use is by Oh and Kim (O-K Model) [ij , 
which is an isobar model. This model starts with Born 
terms, which include t-channel (with K, K* and k ex- 
changes), s-channel ground state nucleon exchanges and 



it-channel A, E and E* exchanges. Additional s-channel 
nucleon resonance exchanges were added to the model 
using the known resonances from the PDG in Ref. [20| . 
referred to here as the K-N-O-K model. One attractive 
point of these models is the inclusion of diagrams with a 
light k meson exchange in the t-channel. As mentioned in 
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TABLE II: Summary of systematic uncertainties. 



K*~ 


"A channel K*^ 


~V° rViannpl 


Event Selections 


2.9% 


4.5% 


Peak Fitting 


7.4% 


5.8% 


Detector Acceptance 


9.2% 


5.7% 


Beam Flux 


7.0% 


7.0% 


Other Sources 


2.5% 


2.5% 


Total 


14% 


12% 



the introduction, the k meson has not yet been firmly es- 
tablished, and these models allows us to study the effect 
of possible k exchange. 

The other model shown here is the Ozaki, Nagahiro 
and Hosaka (O-N-H) Model [2l[, which is a reggeized 
model. This model takes into account all possible hadron 
exchanges with the same quantum numbers (except for 
the spin). The coupling constants and k exchange pa- 
rameters are the same as those used in O-K Model [4|. 

Fig.HOlshows those calculations compared with our dif- 
ferential cross sections, where the solid curves represent 
the theoretical calculations from the K-N-O-K Model and 
the dashed curves represent the O-N-H Model. The cor- 
responding curves are shown in Fig. [TT] for the total cross 
sections, where the curves are explaind in the figure cap- 
tion. The O-K model includes s-channel diagrams with 
most well-established nucleon resonances below 2 GeV 
whereas the K-N-O-K model includes two additional 
s-channel resonances up to 2.2 GeV. Interpretation of 
these results are discussed in section IVTl 

Fig. shows the total cross section ratio of the re- 
actions -yp -> K*°Y>+ to 7p ->■ K*+A. The K*+A data 
alone are not very sensitive to the n exchange due to 
the unknown strength of the coupling constant, g K NA- 
However, the coupling constants of these two reactions 
is related in the effective Lagrangian models, and so the 
ratio is sensitive to the effects of k exchange. The dots 
with error bars in Fig. [T^] use the present data along with 
the previously published CLAS data for K*°E+ @. We 
note that another data set exists for the K*°Y> + reaction 
from CBELSA[6|, but we have chosen to use CLAS data 
in both numerator and denominator to reduce systemat- 
ics. The two curves are the theoretical predictions from 
O-K Model I and II Q, where Model I includes minimal 
t-channel k exchange, while Model II has a significant 
contribution from re exchange. 



VI. DISCUSSION AND CONCLUSIONS 

We presented here the first high-statistics measure- 
ment of the reactions jp — > K* + A and jp — > K* + T,°. 



The data are from the glla experiment using the CLAS 
detector at Thomas Jefferson National Accelerator Fa- 
cility. Differential cross sections are presented for nine 
equal-spaced bins in cos 0^2"^ for each photon energy bin 
of 0.1 GeV width from threshold (1.7 or 1.8 GeV, re- 
spectively) up to 3.9 GeV. Total cross sections, based on 
fits to the differential cross sections are also presented for 
both reactions. 

The cross sections for the K* + A final state are com- 
pared with calculations from two effective Lagrangian 
models, one based on an isobar model and the other 
based on the Regge model. Neither calculation matches 
the data over the broad kinematic range measured here, 
but the isobar model compares more favorably, especially 
at higher photon energies. However, both models signif- 
icantly underpredict the total cross sections in the range 
2.1 < E 1 < 3.1 GeV. Inclusion of all well-known nucleon 
resonances improves agreement with the data in the re- 
gion of E 1 ~ 2 GeV, but has only a small contribution 
above ~ 2.3 GeV, and cannot explain this excess cross 
sections in the new data. 

It remains an open question whether the excess 
strength of the K* + A final state in this photon energy 
region is due to additional couplings to yet-unidentified 
nucleon resonances at higher mass, or whether it is due to 
other effects such as channel-coupling through final-state 
interactions or interference at the amplitude level with 
other physics processes such as photoproduction of the 
K°T,*+ final state. The latter effect was studied using 
a simplified Monte Carlo generator and showed little or 
no effect due to interference with the K°^* + final state, 
but more sophisticated theoretical calculations should be 
done to study interference effects. 

In comparison, the K* + Y,° final state has a sharper 
peak in the total cross section at W ~ 2.25 GeV, and 
falls off more quickly with increasing photon energy than 
for the K* + A final state. This suggests whatever mech- 
anism that causes the excess cross section for the latter 
final state is not present in the K* + Y.° photoproduction. 
However, theoretical calculations are not yet available for 
this final state, and we must wait for more theoretical de- 
velopment before any such conclusion can be reached. 

One of the goals of this measurement was to under- 
stand the role of the k meson exchange, which can con- 
tribute to K* + photoproduction but not to K + pho- 
toproduction. Although no definite conclusion can be 
reached from the present data, the ratio of total cross 
sections for the K* + A and the K*°T, + final state com- 
pared with a similar ratio calculated in the model of Oh 
and Kim suggests that the model with significant k ex- 
change is in better agreement with the data ratio. This 
agrees with the conclusion from a recent study of the 
beam asymmetry measurement Q of the K*°T, + final 
state using a linearly polarized photon beam at forward 
angles. However, we must be careful in making any firm 
conclusion regarding the role of the k exchange until the 
theoretical models have better agreement with the K * + A 
total cross sections above ~ 2.1 GeV. The excess strength 
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FIG. 10: Differential cross sections of jp — > K* + A plotted for incident photon energies from 1.7 to 3.9 GeV. The solid (blue) 
curves represent the theoretical calculations from the K-N-O-K Model [201 ] and the dashed (magenta) curves represents the 
calculations from O-N-H Model with resonance terms 12111. 



of the new data above 2.1 GeV may change the effects 
of k exchange in the ratio. However, the general idea of 
comparing the K* + A and K*°Yj + cross sections, which 
are affected differently by k exchange, is something that 
can be studied now that these new data are available. 
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FIG. 11: Total cross sections of the reaction 7p —t K* + A. 
The solid (cyan) and dash-dotted (blue) curves represent the 
theoretical calculations from the O-K and K-N-O-K mod- 
els, respectively. The dotted (magenta) and dashed (green) 
curves represents the O-N-H model with and without reso- 
nance terms. 
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K*°T, + to 7p — ► K* + A. The ratio uses the present data 
in the denominator and data from Ref. [f| in the numerator. 
The dashed and solid curves are theoretical calculations from 
Oh and Kim 4] Model I and Model II, respectively. 
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